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ABSTRACT 


The purpose of this paper is to examine long wave propagation 
in a shallow channel of triangular cross section. Solutions for 
small scale (f = 0), large scale (f= 0), symmetric and asymmetric 
channels are obtained. Results are shown to be consistent with 
earlier work for the fundamental mode (M = O), and with edgewave 
solutions over a gently sloping bottom for the higher modes. A second 
class of waves (quasigeostrophic waves) is also obtained when the 
Coriolis effect is included. The results are compared with those for 


a rectangular channel. 
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maximum depth of channel 
gravitational acceleration 

velocity components 

X-variation of velocity components 
surface height 

X-variation of surface height 

slope 

slope factor (ne 1) 

frequency 

Coriolis parameter 

wave number 

transformed lateral distance 
non-dimensional wave number 
non-dimensional frequency (squared) 
Laguerre polynomial with argument a 
arbitrary constants 

wave celerity 

non-dimensional frequency 

wave mode 


channel scale parameter 
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I, INTRODUCTION 


Investigations of long wave propagation in channels with 
cross sections which don't vary along the channel began with 
Kelland [1839], who found solutions for progressive waves ina 
uniform channel of triangular cross section with sides inclined at 
45° to the horizontal. MacDonald f1894] examined wave propaga- 
tion in a triangular channel with sides inclined at 30° to the hori- 
zontal. Most investigators have limited their work to small scale 
channels (i.e. the Coriolis effect was neglected as in the above 
cases) in which velocities normal to the axis of the channel were 
zero. An exception is Dronkers (1964), who studied waves ina 
rotating rectangular channel, allowing cross-channel motion. 

Many rivers, estuaries, and shallow seas can be approximately 
represented by a shallow triangular channel (at least this is more 
realistic than a rectangular channel). The shallow-water solutions 
for long waves progressing along such a channel are closely related 
to the well known solutions for edgewaves [Ursell, 1962; and eid. 
1958 | . oucha solution is carried out in this paper and the results 
compared to the work of Dronkers. Most of the results are obtained 


numerically. When possible, they are compared to asymptotic work. 


II. ANALYSIS 


Long waves are studied, propagating along a channel with a 


uniform triangular cross section. The coordinate system is shown in 


PHeruressl.. 


Oe 7 has 





Figure ] 


h.. The slopes of the channel 


The maximum depth of the channel is 
(under the shallow 


amd ms (nea) for x% 20 


sides are s for x>oO 
The depth then is 


water approximation s,ns<<l ). 
hex) a A nS (x =) (1) 
hx) = ho +nsx (x < 0), (2) 


The equations to be satisfied are the shallow-water equations of motion 


ou of = 

ot = fv ro 5 a 0 (3) 

Ae ake — 

) + + a ees _ =a (4) 
and the continuity equation 
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where u = x -— component of velocity 
V = we- component of velocity 
¢ = surface height 
and y is along the channel. A wave motion in the y direction may 


be prescribed by the substitutions 


f= n(x) ai(my+at) 
ie const? 

; -cmy+ct) (6) 
V = V(x) e 


where mm is the wave number, and o the frequency. Then for o*x a 
_ ig AR 
U =r (mtr + acy, 


V= ap (mons f Fe) 


from equations (3) and (4). Substitution into the continuity equation 


yields 


a 1 
ha) Ga + Gh dn a Ge fh 4 Zat — mih)n = ¢ 


forall x andany h(%). 
Equation (8) will now be dealt with in the regions x >O and x<O 


separately, and the results will be patched to insure wave-height and 


velocity continuity at x =0. 


A. EQUATION FOR X >0 
Let F=a=h./s =x . een d/Ax = - VAE, I 7/y,: = Amd and 


h@) = h.- SX =S$£ . The solution region is shown in Figure 2. 


a! 





SOLUTION 





REGION 
Figure 2 
Then (8) becomes 
d“n *” a 
ian +h 4 (p-m'£)n = 0 (9) 
‘ 2 
where P ioe ae 7 .. F ug 


B. EQUATION ROR’ X < 0 


Let €= X + penis . Then h@&)=h.+nsx=ns§ and the corres- 


ponding solution region is shown in Figure 3. 





SOLUTION 
REGION 





Figure 3 
Then (8) becomes 
' A” AN a 2! - 
Pap * Se + (P m*' )n= 0 (11) 
a 
where p: = ae 4 ant (12) 
g@ns oa . 


uz 


c. SOLUTION 


Using the substitutions 





ee) = a ea Cx > 0) 
n@y- E786 (x90) ” 
equation (9) and (11) are put into hypergeometric form (Slater f1960) ): 
acc 
ac ~W ae La) ” 
dwi + (iW) e = 0 (x 70) (14) 
a ') 
V' a& oe 6 / Pe P71) (15) 
aye + Alas Gi = 0 (x < a) 
where W= 2€™, and w= gm : 
These equations have solutions of the form 
G= AFfaCw)+ GH(w) Cx >0) 
(16) 
= B Fa (w') + DHW') (x <0) 
where 
— PS 
am 
Bee p' -- om eo 
~ cA 
But H, H (U in Slater's notation) are unbounded as § and ¢ 
become small and must be discarded. Then the solutions of (14) and 
(15) become 
~me 
m(f)=Ahe™ * Fa WwW) Cx >0) _ 
: =n e" ‘ 7 
n(@)=Be Fai (Ww) CEO) (19) 


1 


where 


[Pau = | © fay eee a" 1) ae a aes) Cs) apo aa 


al(a'+ 1) 12 gi(a'+i)(a'+ a) 3 


\ ae La 
fa Gv) = 1 + aw + ane + —. ee ae 


The quantities A and B are constants; we can set A = 1 as the 
results will be independent of wave amplitude in this linear theory. 


At xX =0 the patching conditions are: 


n(E)= n(E') at x =0 (W=eWe, = wi) (20) 


(i.e. , the wave-height is continuous), and 


U = U’ at x=0 (We We, W'= We ) (21) 


(i,.e., the cross-channel velocity component is continuous). Note 
that We = amh./s and \y,’ = am h./ns = We/n are convenient 
nondimensional wave numbers. 


Using (20), «2 1)becomes 


AN - AN 
Se -_ s erga (22) 
or 
(Et (Ble 
Alf Behe Ag’ go hh. 0. (23) 


Then, using (7), (20), and (23), V is also continuous at x= 6 
For fa (we) = 0 , (20) gives: 1) B= O- which means (21) 
Cannot be sarictied Or 2m FG) 50 which corresponds to a 
node at X=0O. For FF. (W.) x O, equations (20) and (23) give, 


eliminating §&, 
14 


i 7 \ Ate 
F. (we) Fate) ~ Fa (we) $28 (be) - Far (A) SE (we) = 0 (2a) 
Note that the solutions f. (Ww.)= Fa’ (w.') = O are contained in 


(24). Also, when the channel is symmetric (n=1.0), (24) becomes 


F. (w.)) Fe (Ww) - 2 ae (w.)] = fe). (25) 


1S 


TI. SMALL-SCALE CHANNELS 


For relatively short waves in narrow channels, where the 


Coriolis effect may be neglected, (10), 112), and @7) become 











and 
o> 
wa s 07 P— da 
Equation (27) can be written 
oth, _ wy, 
x = gst = (1 1a) 


\ 
where of is a convenient nondimensional frequency. Similarly, 


a 


a ’ 
oS Ga ' aa 
2h 
SC A Ste = 1) We 
and “Hh sia gns* = (1 2a) a] ° 


Equating like terms in (28) and (30) yields the following relationship 


between a and a’: 
(1- 2a) = n(i- Qa’). 


Then equation (24) relates nondimensional wave numbers \W, 


frequencies &. 


A. SYMMETRICAL CHANNEL (f= /.0) 


i ig Asymptotic Results 


For channels with sides of equal slope (31) gives a 


If we<¢<< 1 =, and a is bounded, an asymptotic result can be 
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(26) 


(28) 


; 
a. 


obtained. Dropping terms of OiGah” ) and higher in (25) gives 





(l-JQa)= WA | Jas a(a~')] (ioc 2%) (32) 
so that ax 4 . Thus, using (28), (32) gives 
we = amh. = Yo * ( é 1) 
GS ~ wags aS (33) 
or 
= 
Since wave celerity ¢ = s/m , this can be written 
Ya qh (35) 


which is identical to Kelland's result for a symmetric triangular channel. 

Since the average depth is ho /2. , then this result also corresponds 

to the classical result for long waves over a flat bottom: c*= q h . 
Other asymptotic results for \4Z<<1 and a unbounded can 

be obtained by iteration. An example is given below. 


Let Ww be the root obtained by keeping terms through \/” 


in Fa (Ww) = 0 Then: 
-aw, = | 
. guy, = | # ae Ww ~ l + ~onF + O(4) 
Lawns | + ate + SGA) wy,’ 
= 1 + phi + age] cis fi + chp] 0(2) 
—- AW = ete 
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This converges rapidly to 
4aWwe=- “Grado 


and from (28), 


o{ «Gzeen/enbee . (36) 


Similar techniques could give the other asymptotes suggested by the 
numerical results presented below. 

For W.>>1 (short waves), it is expected on a physical basis 
that the waves in the two parts of the channel become independent of 


one another, and that Ursell's edgewave result 


ay eo GM hyo.) (37) 


would hold. Although this has not been confirmed here analytically, 
it is borne out by the numerical results. 

2% Numerical Results 

A Fortran IV computer program to solve (24) for small scale 

channels (f = o) is given in Appendix A. The results for the 
symmetric channel (mM =/.0) are shown in Figures 4 and 5. In Figure 5, 
the dispersion relation for the fundamental mode (M = 0) and the 
next five modes are plotted. 

For \/4 << 4 (Figure 4) the calculations agree well with the 
asymptotic results (34) and (36). 

For , Wee > >sel: , note that 


xo = (aM +i) CM 


) 

oS OO 
9 YW 
JY nL 


or me x Camel) gons CM 
which is the same as (37). 
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Dispersion Relation: Symmetric Channel (W << 1 ) 
(f=0,n=1.0) 


as, 


100.0 
. NS ae 
10.0 


$ 1.0 £.0 10.0 50.0 
We 
Figure 5 
Dispersion Relation: Symmetric Channel (f=0,n=1.0) 
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3m Wave Profiles 
Figure 6 is a sketch of the wave profiles for modes 

0 through 5 showing the nodal points for the symmetrical channel. 
Note that the nodes for all modes are symmetric about the center of 
the channel. The profile for the fundamental mode (M= 0) has no 
nodes along the channel. Wave mode (M=1) is the first asym- 
metric profile, corresponding to the solution n(£) =— 7 (8'), [Fa (we ) = o}. 
Wave mode (M =2) has a symmetric profile with two nodal points. 
The higher modes have symmetric profiles for even mode numbers 
and asymmetric profiles for odd (M), with the number of nodal points 
equal to (M). The nodal points (except x =o ) move toward the 
channel walls with increasing \J,. This corresponds to the edge- 


wave condition for short waves (large ™ ). 


B. ASYMMETRICAL CHANNEL (rn *& /.0) 


d. Asymptotic Results 


For \4<< i and a bounded, (24) can be written 


a coc) = Ww.[aa’ (Dit) +: ana) 4. — 0] a 0 (w:) ! 





(38) 
Since W,<<1, then atazxi. Using (31), this gives 
a+ zp, (aat n-1)= | 
so that 
Gye; al ae 
Let a=Yy -€ and al= i ae Then (38) gives 
orn t 
w= poo (e+e) (39) 





=e 


Figure 6 


Wave Profiles: Symmetric Channel (n=1.0) 


which with (28) and (30) gives 


2 


ie = ete, 


This is the same as (34), and is independent of Nn. 
Other results for \/, <<]. similar to (36) are expected for 
asymmetric channels but are not examined here. 
Ze Numerical Results 
The results for n= 2.0 and mn =5.0 are shown in 
Figures 7 and 8. In each case, the asymptotic result (34) is confirmed. 
In fact, the entire dispersion relation for the fundamental mode (M= 0) 
appears to be independent of Nn. The higher modes have all been 
shifted to higher frequencies in comparison with the symmetric case. 
ce Wave Profiles 
Profiles for n=2.0 and n= 5.0 are shown in Figures 9 
and 10. Again the fundamental mode (M=0) is unchanged, but the 


higher modes have been altered significantly from the symmetric case. 


Zs 


Dispersion Relation: 


5.0 {0.0 50.0 
Ve 
Figure 7 


Asymmetrical Channel (f=0,n=2.0) 


10a0 


$0.0 “ 
120 
5.0 
4 
4.0 
.§ 
ef 
5 4.0 


5.0 10.0 50.0 
Whe 


Figure 8 


Dispersion Relation: Asymmetrical Channel (f=0,n=5.0) 


IAS 





Figure 9 


Wave Profiles: Asymmetric Channel (n=2.0) 





Figure 10 


Wave Profiles: Asymmetric Channel (n=5.0) 


fay 


IV. LARGE-SCALE CHANNELS 


For large channels where the Coriolis effect must be considered, 


(10) -(12)) and (17).gige 
_ Das ae 


i 
jf = fa = gms 1) — op 
and (40) 


(wl) + ob 


a SS 20’ = of 


where = o/f ‘ 


Equating like terms in (40), 


ie oe (41) 
and 
= Wile I BES): be _ (42) 
where K = af he 


Now W must satisfy both (24) and (42). 


A. SYMMETRICAL CHANNEL (n= /.0) 
1s Asymptotic Results 
From (41), a and a’ arerelated by aq=a’+ VANe But 
for We. <<1 and uw bounded, (42) implies that Ja +a'| is Targey 


sothat az aq’. Then (40) gives 





¢* cS 
-Ja = (or ry - + (43) 


This can be written 


| 
s 
q 
}e2 
<7 
| 
~ 
is 
Ls 


~la = 


28 


From (36) and (44), 








aH _ 2 oe ea 
ef = (ur? 1) Ae Cm << R) 
Dropping the small term gives 
ae 
a 5 Ye gy (45) 
2. Numerical Results 


A Fortran IV computer program to solve (24) and (42) is given 
in Appendix B. A value of , =0.! is chosen, which corresponds approxi- 
mately to the dimensions of the Red Sea. The results for the symmetric 
channel are shown in Figure 11. Note that for \A << 1], the asymp- 
totic result (45) is confirmed. For Ur << 1{ anewclass of waves 
appears. These are analogous to the quasigeostrophic edgewaves 
obtained by Reid. Figure 12 is an expanded plot of these low frequency 
waves. 

For GW, We.>>1 #£(short high-frequency waves), the Coriolis 
effect should be negligible. Then, fora fixed WA, the frequencies o 


for f=0 and f*°0should be the same. Thus, 


od SES h. o-* 
3. eh/e ; 


a 


This relation is observed to hold for wyW.>>1. 


B. ASYMMETRICAL CHANNEL (N23 2.0) 


ies Asymptotic Results 


For W<<1 and a _ unbounded, asymptotic results similar 


to (45) are expected, but are not derived here. 


a 








peseeeeeeiereee 
: a r 6 > 10 
\ws 
Figure i11 
Dispersion Relation: Symmetric Channel (f*0,n=1.0) 
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La Numerical Results 





The results for n=2.0 and nms=5.0 are shown in Figures 
13 and 14. Note that all modes have been shifted to higher frequencies 
in comparison with the symmetric case. Also for UF, \4>>@ We 


f=0 and fxo results give good agreement via (46). 


a2 





+0 


Dispersion Relation: 


6 g 10 
We 
Figure 13 


Assymetric Channel (f0,n=2.0) 
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4 
a) 
i 2 4 “6 & 10 
\We 
Figure 14 
Dispersion Relation: Asymmetric Channel (fx0,n=5.0) 
34 


V. COMPARISON OF RESULTS FOR A RECTANGULAR CHANNEL 
AND A SYMMETRIC TRIANGULAR CHANNEL 


The dispersion relation due to Dronkers for a channel of uniform 


depth (h) and width (Ww) is (Appendix C) 


a a 
a - 41 Co Oren.) =~ 


It is reasonable to compare channels of equal average depth, (Figure 15). 


\/ = ah. S$ 






SLOPE =-S—”? Nest oPEe =5s 


Figure 15 


Then h=h./a and W= 2h./'s , and (47) becomes 





W. = a (7 *~ £7) - M*rr* oe) 


~e SMALL-SCALE (f = 0) 


For small-scale channels (f=0), (48) becomes 


Wee eme nn §=—S (0, 1, 2 (49) 


This agrees with (34) for the fundamental mode (M =0). This dispersion 
relation (49) is compared with that for a symmetric triangular channel 


in Figure lo. 


So 


'6a.Q 


§a0 


10.0 


5.0 


Rectangular Channel 


Triangular Channel 





5 Le 5.0 16.0 $6.0 
VW, 
Figure 16 
Dispersion Relations: Symmetric Triangular 


vs. Rectangular (f=0) 
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For large-scale channels (f ¥ 0), (48) becomes 
a 
Wer a4 KK (wi 1 -M rr (mM = Ouleuccn ey (S50) 


This dispersion relation (50) is compared with that of a symmetric 


triangular channel in Figure 17. 
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Dispersion Relations: 





Rectangular Channel 


Triangular Channel 


6 $ 


We 
Figure 17 


Symmetric Triangular 


vs. Rectangular (fx0) 
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APPENDIX A 


FORTRAN IV PROGRAM TO SOLVE PATCHING EQUATION 
FOR SMALL-SCALE CHANNELS 


IMPLICIT REAL*8 (A-H , O-Z) 
DIMENSION C (36) ,C (36), CP(36) , DP(36) ,CCP(36) ,CDP(36), DC P(36), WORK 
* (36) ,B(36) , ROOTR(36) , ROOTI (36) 
READ(S,5)A 

FORMAT (F10.5) 
IF(A.EQ.-9.50)GOTO99 

S—5), 0 

AP=(S+2.0*A-1.0)/(2.0*S) 
WRITE(6,10)A,AP,S 

meme ros ,'A=',F10.5,5X,'AP=', F10.5,5X,'N=', Pome) 
C(1)=1.0 

Giay=A 

WO 2 N =3,12 

C (N)=C (N-1)* (A+ (N-2)) /((N-1)* (N-1)) 
CONTINUE 

Pip =A 

D(2)=A* (at+1)/2.0 

me 32 N=3,12 

D(N)=D(N-1)* (A+(N-1))/((N-1) *N) 
CONTINUE 

GR Se 

CP(2)=AP/S 

WO 4 N=3,12 

CP(N)=CP(N-1)* (AP+ (N-2)) /(N-1) * (N-1)*S) 
CONTINUE 

DP(1)=AP 

DP (2)=AP* (AP+1)/(2.0*S) 

Were Nac, 12 

DP(N)=DP(N-1)* (AP+ (N-1))/((N-1)*N*S) 
CONTINUE 

m@7 K=112 

Gerij=0.0 

ee 1—1, 12 

Wee, j= 12 

IF(I+J.NE.K)GOTO7 

CCP (K-1)=CCP(K-1)+C (I) *CP (J) 
CONTINUE 

De ca), 12 

GDrig=0.0 

Dero. 12 

Bers jJ=l,12 

IF(I+J.NE.K)GOTO8 
CDP(K-1)=CDP(K-1)+C (I) *DP(J) 

39 


Tal 


Zo 


30 


99 


CONTINUE 
YO) Keay 

DEPiIG—0. 0 

DO 9 1=leae 

DOs) 12 

IF(I+J. NE.K)GOTO9 

DC P(K-1)=DC P(K-1)+D (I) *C P(J) 

CONTINUE 

DO Me es) 

B(N)=CCP(N)=CDP(N) -DCP(N) 

CONTINUE 

CALL DPOLRT(B, WORK, 10, ROOTR, ROOTI, IER) 
WRITE (6,25) (ROOTR(I) , ROOTI(I) , I=1, 10) 
FORMAT (//,4X,2D25.16,/ 
WRITE(6,30)IER 

FORMAT (50X, I2) 

GOTO] 

STOP 

END 


40 
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APPENDIX B 


FORTRAN IV PROGRAM TO SOLVE PATCHING EQUATION 
FOR LARGE-SCALE CHANNELS 


IMPLICIT REAL*8(A-H , O-2) 


DIMENSION C (36) ,D(36) ,CP(36) , DP(36), CC P(36), CDP(36), DC P(36), WGRK 
#90) ,B(36) ,ROODTR(36) , ROOTI(36) 


er—-1.0p=t2 
H=1000.0 
G=9.8 
SLP=0.1 

5=5 .0 
CAPA=.10 
READ (5 ,5)W 


FORMAT (F10.5) 

IF(W.GT.14.5)GOTO99 

WRITE(6,10)W 
FORMET (Sen, V=",FIO.5 ,/) 
OMEGA=C .50 

OMEGA=OMEGA-.001 
IF(OMEGA.DE.0.0)GOTO!1 
A=.5*(1.0/OMEGA+1. 0-(CAPA* (OMEGA* OMEGA-1.0)) /W) 
AP=.5*(1.0-1.0/OMEGA-(CAPA* (OMEGA* OMEGA-1.0)/(W*S))) 
Cig)=1.0 

Go =A 

me 2 N=3, 12 

C (N)=C (N-1)* (A+(N-2)) Z((N-1)* (N-1)) 
CONTINUE 

Dilj=A 

D(2)=A* (A+1)/2.0 

we seN=3,12 
D(N)=D(N-1)* (A+ (N=-1))/((N-1) *N) 
CONTINUE 

Chi =1,0 

CP(2)=AP/S 

we 4 N=3, 12 

CP(N)=CP(N-1)* (AP+(N-2))/((N-1)* (N-1)*S) 
CONTINUE 

maj =AP 

DP(2)=AP* (AP+1) /(2.0*S) 

wo 6 N=3 12 

DP(N)=DP(N-1)* (AP+(N-1))/((N-1)*N*S) 
CONTINUE 

De. 7 K=1,12 

Serih)=0.6 

DO 7 I=1,12 

me 7 j=1 12 

IF (I+J. NE.K)GOTO7 4] 


CCP(K-1)=CC P(K-1)+C (I) *CP(J) 
7 CONTINUE 
DO 8 Keil, 12 
CDP(K)=0.0 
DOPaM= 1,12 
DO 8 J, 12 
IF(I+J. NE. K)GOTO8 
CDP(K-1)=CDP(K-1)+C (I) *DP(J) 
See @ONMNLIE 
Derg 1-12 
DCP(WEHC.0 
WOro j= 12 
PO jal 12 
IF(I+J).NE.K)GOTO9 
DC P(K-1)=DC P(K-1)+D (I) *CP(J) 
9 CONTINUE 
m@e1)] N=1,11 
B(N)=CC P(N) -CDP(N) -DCP(N) 
11 CONTINUE 
FC N=B(1)+B(2) *W+B (3) * W* W+B (4) * W**34+B (5) *W**4+B (6) *W**54B(7) *W** 
*+B(8)*W**7+4B(9)*W**8+B(10)*W**94+B(11)*W**10 
WRITE(6,15)A,AP, OMEGA, W,FCN 
15 FORMAT(2X,4F15.7,D25. 16) 
GOTO12 
99 STOP 
END 
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APPENDIX C 


DISPERSION RELATION FOR DRONKERS' RECTANGULAR CHANNEL 


A. FORMULATION 


The equations to be satisfied are: 


Substitution of 


(u,v, if) = (u, v, ») rapa le 


into (1) and (2) gives 
icU-fv4 ggt=o 
aay, a tera gre O 


n(g + imv) econo, 





Then, for o*x f* , 
een, Ax 


Ve aai(meont f Am 


and the ordinary differential equation for 7 is 





An 2 
dx —° 
where re Z Sy a ; ae 


with solution of the form 


N= a aimAx+ brow Ax, 
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(1) 


(2) 


(4) 


Be BOUNDARY CONDITIONS 


The coordinate system is shown in Figure Cl. The flux normal 


to the channel walls must vanish at the walls. 


xX =\/ EY AO A ME ME SO LE ML Me Se Lee ey ed Se Ow 





Figure Cl 


Then (4) gives 


min 4 dm _ g 


z (x= OW) . 


Substitution of (7) into (8) yields 
ed AW 
| ae + BY Item AW = 0 


so that, 


tan AW = 0 


or Ar = Ae Gl =n), 


Then (6) gives the dispersion relation 


2 a om 2. 
m* = me il fr 
qh al 
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(11) 
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